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Thermal axion production in the primordial quark-gluon plasma 

Peter Graf and Frank Daniel Steffen 
Max- Planck- Institut fur Physik, Fohringer Ring 6, D-80805 Munich, Germany 

We calculate the rate for thermal production of axions via scattering of quarks and gluons in the 
primordial quark-gluon plasma. To obtain a finite result in a gauge-invariant way that is consistent 
to leading order in the strong gauge coupling, we use systematic field theoretical methods such 
as hard thermal loop resummation and the Braaten-Yuan prescription. The thermally produced 
yield, the decoupling temperature, and the density parameter are computed for axions with a mass 
below 10 meV. In this regime, with a Peccei-Quinn scale above 6 x 10 8 GeV, the associated axion 
population can still be relativistic today and can coexist with the axion cold dark matter condensate. 

PACS numbers: 14.80.Va, 98.80.Cq, 95.35.+d, 95.30.Cq 
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I. INTRODUCTION 

If the Peccei-Quinn (PQ) mechanism is the explana- 
tion of the strong CP problem, axions will pervade the 
Universe as an extremely weakly interacting light parti- 
cle species. In fact, an axion condensate is still one of the 
most compelling explanations of the cold dark matter in 
our Universe [l], • While such a condensate would form 
at temperatures T < 1 GeV, additional populations of 
axions can originate from processes at much higher tem- 
peratures. Even if the reheating temperature Tr after 
inflation is such that axions were never in thermal equilib- 
rium with the primordial plasma, they can be produced 
efficiently via scattering of quarks and gluons. Here we 
calculate for the first time the associated thermal produc- 
tion rate consistent to leading order in the strong gauge 
coupling g s in a gauge-invariant way. The result allows 
us to compute the associated relic abundance and to esti- 
mate the critical Tr value below which our considerations 
are relevant. For a higher value of Tr, one will face the 
case in which axions were in thermal equilibrium with the 
primordial plasma before decoupling at the temperature 
Td as a hot thermal relic [l], 0] . The obtained critical Tr 
value can then be identified with the axion decoupling 
temperature To- 

We focus on the model-independent axion (a) interac- 
tions with gluons given by the LagrangiarQ 



£a = 



9s j^ib Pyb ii v 

32tt2/pq ^ ' 



(1) 



with the gluon field strength tensor G v , its dual G b = 
£ l iupoG bf " J /2, and the scale /pq at which the PQ sym- 
metry is broken spontaneously. Numerous laboratory, 
astrophysical, and cosmological studies point to 



/ PQ >6x 10 8 GeV 



(2) 



1 The relation of /pq to the vacuum expectation value (<j>) that 
breaks the U(1)pq symmetry depends on the axion model and 
the associated domain wall number N: /pq oc (<f>)/N; cf. U01 
and references therein. 



which implies that axions are stable on cosmological 
timescales 0, Q • Considering this /pq range, we can ne- 
glect axion production via 7T7t — > ira in the primordial hot 
hadronic gas @, 0j- Moreover, Primakoff processes such 
as e~7 — > e~a or 177 — > q a are not taken into account. 
These processes depend on the axion model and involve 
the electromagnetic coupling instead of the strong one 
such that they are usually far less efficient in the early 
Universe 

We assume a standard thermal history and refer to 
Tr as the initial temperature of the radiation-dominated 
epoch. While inflation models can point to Tr well above 
10 10 GeV, we focus on the case Tr < /pq such that no 
PQ symmetry restoration takes place after inflation. 

Related studies exist. The decoupling of axions out 
of thermal equilibrium with the primordial quark-gluon 
plasma (QGP) was considered in Refs. [3, 0|. While the 
same QCD processes are relevant, our study treats the 
thermal production of axions that were never in thermal 
equilibrium. Moreover, we use hard thermal loop (HTL) 
resummation [§] and the Braaten-Yuan prescription [l(| , 
which allow for a systematic gauge-invariant treatment 
of screening effects in the QGP. In fact, that prescription 
was introduced on the example of axion production in a 
hot QED plasma [l(J; see also Ref. [llj ]. 



II. THERMAL PRODUCTION RATE 

Let us calculate the thermal production rate of axions 
with energies E > T in the hot QGP. The relevant 2^2 
scattering processes involving (JlJ are shown in Fig. [T] 
The corresponding squared matrix elements are listed in 
Table H where s = (Pi + P 2 ) 2 and t = (Pi - P 3 ) 2 with 
Pi, P2, P3, and P referring to the particles in the given 
order. Working in the limit T ^> nij, the masses of all 
particles involved have been neglected. Sums over initial 
and final spins have been performed. For quarks, the 
contribution of a single chirality is given. The results 
obtained for processes A and C point to potential IR di- 
vergences associated with the exchange of soft (massless) 
gluons in the t channel and u channel. Here screening 
effects of the plasma become relevant. To account for 



2 



Process A: g a + g b — >• g c + a 



9" 



+ 





1j / V) 9" 

Process C: qi + g a — !> qj + a (crossing of B) 

FIG. 1. The 2 — >■ 2 processes for axion production in the QGP. 
Process C also exists with antiquarks q~ij replacing qtj. 



such effects, the QCD Debye mass tod = y3m g , with 
m g = g s TyjN c + (rt//2)/3 for N c = 3 colors and u/ = 6 
flavors, was used in Ref. Q. In contrast, our calculation 
relies on HTL resummation [t| [13] which treats screening 
effects more systematically. 

Following Ref. (loj . we introduce a momentum scale 
fe cu t such that g s T <C fc cut <C T in the weak coupling 
limit <? s <C 1. This separates soft gluons with momentum 
transfer of order g s T from hard gluons with momentum 
transfer of order T. By summing the respective soft and 



hard contributions, the finite rate for thermal production 
of axions with E > T is obtained in leading order in g s , 



E dWa^ E dW a 



d 3 p 



d 3 p 



soft + E d3 P 



(3) 



hard 



which is independent of fc cut ; cf. ([5} and j7| given below. 

In the region with k < k cut , we obtain the soft con- 
tribution from the imaginary part of the thermal axion 
self-energy with the ultraviolet cutoff fc C ut, 



E 



dW a 



d 3 p 
= Ef B (E) 



1)T 



8192^s/2 Q 



In 



^cut 



fe<fe C ut 



1.379 



(4) 
(5) 



with the equilibrium phase space density for bosons 
(fermions) fB(F){E) = [exp(E/T) =F Our derivation 
of ||SJ| follows Ref. [l(| • The leading order contribution to 
ImII a for k < fc cut and E >T comes from the Feynman 
diagram shown in Fig. O Because of E > T, only one 
of the two gluons can have a soft momentum. Thus only 
one effective HTL-resummed gluon propagator is needed. 

In the region with k > fc cu t, bare gluon propagators 
can be used since fc cut provides an IR cutoff. From the 
results given in Table |T] weighted with appropriate mul- 
tiplicities, statistical factors, and phase space densities, 
we then obtain the (angle-averaged) hard contribution 



E 



,dW„ 



d 3 p 



1 



dCl r , 



d 3 p 3 



n (2703225, 



hard 2(27T)3./ 4tt 

x (2n) i 5 4 (P 1 +P2-P3- P)Q(k - fccut) 
x ]T /i(£i)/ 2 (£ 2 )[l ± / 3 (^ 3 )]|M 1+2 ^ 3+a | 

9^ - 1) 



= E- 



5127T 7 /| 



PQ 



/ B (i5)r 3 
n f j^z ln ( 2 ) 



48tt 



c 2 J 48tt 



1)1 ( p- ) - ' ^ ' 



17 

t 



(6) 



2C'(2) 



+^c(4bB 



7-< 3 ) ^ j-t. C7-W , r( 3 ) A 
j bbbJ "t" "A-'fbf "t j ffb/ J 



C(2) 
(7) 



TABLE I. Squared matrix elements for axion (a) production 
in 2-body processes involving quarks of a single chirality (qi) 
and gluons (g a ) in the high-temperature limit, T>mi, with 
the SU(A r c ) color matrices / ai>c and T 3 ". Sums over initial and 
final state spins have been performed. 

Labeli Process i \Mi\ 2 / ( t^zjt^ J 

A g*+g>^g°+ a -4 *^±g±g£ |/° 

B gi + fc^ fl a +a (2£+2t + s ) |T/ t | 2 

C * + 0°-»-fc+a (_3d_ 28 -t) |7-| 2 



.9 




FIG. 2. Leading contribution to the axion self-energy for soft 
gluon momentum transfer and hard axion energy. The blob on 
the gluon line denotes the HTL-resummed gluon propagator. 
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with Euler's constant 7, Riemann's zeta function £(z), 



1 



'bbb(fbf) 327r 3 



dEo 



x -G(Ei - E 3 




d 
dE[ 



E+E 3 



dE x In 



E 2 I rp 2 



\E\ — E3 

e7~ 



/bbb(fbf)-^-(£'i + E 3 ) 



+ Q{E 3 - £'l)^j-[/BBB(FBF)(£ ; i + El)} 



1 



2 

2f2 



7bbb(fbf)I — g2 3 J J 



(8) 



rO) 

'bbb(ffb) - 327r 3 
x f 0(£ - E 3 



dE 3 



e+e 3 



dE 2 f } 



BBB(FFB) 



EM 



Q(E 3 -E)- 



E 2 (E 3 +E) ' v J 'E 3 + E 
[Q(E 3 - E)Q(E 2 - E 3 ) - 0{E - E 3 )Q{E 3 - E 2 )\ 
E2 — E 3 



E 2 



[E 2 (E 3 - E) - E 3 (E 3 + E)} 



/BBB,FBF,FFB — fl(Ei)f2\E2)[l ± f3\E 3 )]- 



(9) 
(10) 



The sum in ([5]) is over all axion production processes 
1 + 2 — > 3 + a viable with ([1]). The colored particles 1-3 
were in thermal equilibrium at the relevant times. Per- 
forming the calculation in the rest frame of the plasma, 
fi are thus described by /f/b depending on the respec- 
tive spins. Shorthand notation (ITU1) indicates the corre- 
sponding combinations, where + (— ) accounts for Bose 
enhancement (Pauli blocking) when particle 3 is a bo- 
son (fermion) . With any initial axion population diluted 
away by inflation and T well below T^> so that axions are 
not in thermal equilibrium, the axion phase space density 
f a is negligible in comparison to /f/b- Thereby, axion 
disappearance reactions (oc f a ) are neglected as well as 
the respective Bose enhancement (1 + f a 1). Details 
on the methods applied to obtain our results (JT])-© can 
be found in Refs. [llHll. 



III. RELIC AXION ABUNDANCE 

We now calculate the thermally produced (TP) axion 
yield = n a /s, where n a is the corresponding axion 
number density and s the entropy density. For T suffi- 
ciently below Td, the evolution of the thermally produced 
n a with cosmic time t is governed by the Boltzmann equa- 
tion 



dn a 
~dT 



3Hn a 



3 dW a 



d^p 



d 3 p 



W a . 



(11) 



Here H is the Hubble expansion rate, and the collision 
term is the integrated thermal production rate 



W a 



C(3)ff A 6 r 6 

647T7/2 



PQ 



In 



m 2 , 



+ 0.406 



Assuming conservation of entropy per comoving vol- 
ume element, (fTTj) can be written as dYj v jdt — W a /s. 
Since thermal axion production proceeds basically dur- 
ing the hot radiation-dominated epoch, i.e., well above 
the temperature of radiation-matter equality T ma t =ra d, 
one can change variables from cosmic time t to tempera- 
ture T accordingly. With initial temperature Tr at which 
y a TP (Tri) = 0, the relic axion yield today is given by 



Y, 



TP 



V TP (T 
1 a V J -mat= 



rad ) 



dT 



18.6g 6 s In 



1.501\ /10 10 GeV 



9s 



fpQ 



Wg{T) 

Ts{T)H(T) 



10 10 GeV 



(13) 



This result is shown by the diagonal lines in Fig.[3]for cos- 
mological scenarios with different Tr values ranging from 
10 4 to 10 12 GeV. Here we use g s = g s {Ts.) as described 
by the 1-loop renormalization group evolution [j| 



llN c -2n f /3h 
24tt 2 V Mz 



-1/2 



(14) 



where g 2 (M z )/(4ir) = 0.1172 at M z = 91.188 GeV. 
Since the methods allowing for the gauge-invariant 

derivation require g s <C 1, (|13|) is most reliable for 
Tp, ^> 10 4 GeV. If g s is too large, unphysical negative 
values of ^ are encountered for E < T that lead to an 
artificial suppression of (IT3l) as indicated by the logarith- 
mic factor; cf. Sec. 3 in Ref. [3] where a similar behavior 



10" 1 1 1 1 1 1 1 

Kfo f» 2.6 x 1CT 3 




10 

T R [GeV] 

FIG. 3. The relic axion yield today originating from thermal 
processes in the primordial plasma for cosmological scenarios 
characterized by different Tr values covering the range from 
10 4 to 10 12 GeV. The dotted, dash-dotted, dashed, and solid 



(12) lines are obtained for / PQ = 10 9 , 10 10 , 10 11 , and 10 12 GeV 
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is discussed in more detail. Nevertheless, we think that 
the applied methods still allow presently for the most 
reliable treatment since they respect gauge invariance. 

Note that (fTB")) is only valid when axion disappearance 
processes can be neglected. In scenarios in which Tr ex- 
ceeds Td, this is not justified since there has been an early 
period in which axions were in thermal equilibrium. In 
this period, their production and annihilation proceeded 
at equal rates. Thereafter, they decoupled as hot ther- 
mal relics at Tb, where all standard model particles are 
effectively massless. The present yield of those thermal 
relic axions is then given by Y" a eq = n cq /s s» 2.6 x 10~ 3 . 
In Fig.[3]this value is indicated by the horizontal lines. In 
fact, the thermally produced yield cannot exceed Y° q . In 
scenarios with Tr such that (fT3|) turns out to be close to 
or greater than Y° q , disappearance processes have to be 
taken into account. The resulting axion yield from ther- 
mal processes will then respect Y£ q as the upper limit. 
For example, for /pq = 10 9 GeV, this yield would show a 
dependence on the reheating temperature Tr that is very 
similar to the one shown by the dotted line in Fig. [3] The 
only difference will be a smooth transition instead of the 
kink at y Q T P = y^eq 

IV. AXION DECOUPLING TEMPERATURE 

The kinks in Fig. [3] indicate the critical Tr value which 
separates scenarios with thermal relic axions from those 
in which axions have never been in thermal equilibrium. 
Thus, for a given /pq, this critical Tr value allows us to 
extract an estimate of the axion decoupling temperature 
Td • We find that our numerical results are well described 
by 

/ f \ 2.246 

T D «9.6xl0 6 GeV( J5 ^J (15) 

In a previous study the decoupling of axions that 
were in thermal equilibrium in the QGP was calculated. 
When following [3j but including (TT4"]) . we find that the 
temperature at which the axion yield from thermal pro- 
cesses started to differ by more than 5% from Y a eq agrees 
basically with (TT5|) . The axion interaction rate T equals 
H already at temperatures about a factor 4 below (fT"5)) 
which, however, amounts to a different definition of Td- 

V. AXION DENSITY PARAMETER 

Since thermally produced axions have basically a ther- 
mal spectrum also, we find that the density parameter 
from thermal processes in the primordial plasma can be 
described approximately by 

„ ^(p ai0 )2 +m 2yTP/oq s{Tq) ^ J pc (16) 

with the present average momentum (p a ,o) = 2.701 T a .o 
given by the present axion temperature T a ^ — 0.332 T$ ~ 



0.08 meV, where To ~ 0.235 meV is the present cosmic 
microwave background temperature, h ~ 0.7 is Hubble's 
constant in units of 100 km/Mpc/s and p c /[s(T )h 2 } = 
3.6 eV. A comparison of T a Q with the axion mass 
m a ~ 0.6 meV (10 10 GeV//pQ) shows that this axion 
population is still relativistic today for /pq > 10 11 GeV. 

In Fig.[4]the solid, dashed, and dash-dotted lines show 
nI P/cq h 2 for T R = 10 6 , 10 7 , and 10 8 GeV, respectively. 
In the /pq region to the right (left) of the respective 
kink, in which T R < T D (Tp > T D ) holds, f£ P(oq) /i 2 ap- 
plies, which behaves as oc fpQ (/pq) for m a 3> T a Q and 
as oc /pq (/pq) for m a <C T a fi. The gray dotted curve 
shows Sl cq h 2 for higher Tr with Tr > Td and also in- 
dicates an upper limit on the thermally produced axion 
density. Observation of those axions will be extremely 
challenging. Even Sl cq h 2 stays well below the cold dark 
matter density ricDM^ 2 — 0.1 (gray horizontal bar) and 
also below the photon density fl 7 h 2 ~ 2.5 x 10 -5 (gray 
thin line) [f| in the allowed /pq range ([2]). There, the cur- 
rent hot dark matter limits are also safely respected [l5| . 

In cosmological settings with Tr > Td, axions pro- 
duced nonthermally before axion decoupling (e.g., in in- 
flaton decays) will also be thermalized resulting in Q cq h 2 . 
The axion condensate from the misalignment mechanism, 
however, is not affected — independent of the hierarchy 
between Tr and Td — since thermal axion production in 
the QGP is negligible at T < 1 GeV. Thus, the asso- 




/pq [GeV] 

FIG. 4. The axion density parameter from thermal processes 
for T R = 10 6 GeV (solid), 10 7 GeV (dashed) and 10 s GeV 
(dash-dotted) and the one from the misalignment mechanism 
for 8i = 1, 0.1, and 0.01 (dotted). The density parameters 
for thermal relic axions, photons, and cold dark matter are 
indicated, respectively, by the gray dotted line (Q^h 2 ), the 
gray thin line (f2 7 /i 2 ), and the gray horizontal bar (tlcDwih 2 ). 
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dated density Of IS /i 2 ~ 0.15 (9f(/ PQ /10 12 GeV) 7 / 6 0,0, 

liH ] can coexist with flJ F ^ cq h 2 and is governed by the 
misalignment angle 6i as illustrated by the dotted lines 
in Fig. [4] Thereby, the combination of the axion cold 
dark matter condensate with the axions from thermal 
processes, tt a h 2 = Q^ ls h 2 + nJ P ^ cq h 2 , gives the analog 
of a Lee- Weinberg curve. Taking into account the re- 
lation between /pq and m Q , this is exactly the type of 
curve that can be inferred from Fig. |H Here our cal- 
culation of thermal axion production in the quark-gluon 



plasma allows us to cover, for the first time, also cosmo- 
logical settings with T R < T D . 
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